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ABSTRACT 

The existence of an intergalactic globnlar clnster population in the Coma 
cluster of galaxies has been tested using surface-brightness fluctuations. The 
main result is that the intergalactic globular cluster surface density (iVicc) does 
not correlate with the distance to the center of Coma and hence with the en¬ 
vironment. Furthermore, comparing these results with different Coma mass- 
distribution model predictions, it is suggested that Niqq must in fact be zero 
all over Coma. On the other hand, the results for Aigc and the faint end of 
the galaxy number counts (beyond uir = 23.5) are connected. So Aigc = 0 
settles the slope of this function, which turns out to be 7 = 0.36 ± 0.01 down to 
tur = 26.5. 

The fact that Aigc = 0 all over Coma suggests that globular clusters were 
formed only, or almost only, from protogalactic clouds. None, or perhaps very 
few, could have formed in isolated regions. It also seems inappropriate to ad¬ 
vocate a relationship between intergalactic globular clusters and dark matter 
distributions, although it is true that the relationship could still exist but not 
be strong enough to have been detected. Finally, since our conclusion is that 
intergalactic globular clusters do not exist in Coma, accretion of intergalactic 
globular clusters might not be signihcant in galaxy formation and evolutionary 
processes in the Coma galaxies. 

Subject headings: galaxies:clusters:individual (Coma)—galaxies:star clusters— 
galaxies: formation—galaxies:evolution—galaxies: luminosity function 


1. Introduction 


Intergalactic globular clusters (IGCs) are not bound to any particular galaxy but move 
freely in the potential wells of galaxy clusters. Table 1 will help the reader with the notation 
used in this paper. The existence of IGCs has been considered by a number of authors. The 
hrst discussion about them was by van den Bergh (1956), who measured the distance to the 
Abell No. 4 globular cluster (GC). Obtaining a distance modulus of 20.8, he concluded that it 
is therefore an intergalactic tramp. Van den Bergh (1958) estimated that the number of IGCs 
within the Local Group is one third of the total number of GCs. But the hrst real quantitative 
studies of IGCs in galaxy clusters were done by Forte, Martinez, & Muzzio (1982) and 
Muzzio (1987). 

Muzzio (1987) suggested that galaxies in galaxy clusters suffer several dynamical effects 
owing to encounters with neighboring galaxies and to the action of the general galaxy cluster 
held; as a result, galaxies lose some GCs. Since these processes are very sensitive to the 
distribution of the total mass of the galaxy cluster, the study of IGCs may help us to solve 
the problem of how the missing mass is distributed. In this scenario, the number density 
of lost IGCs should follow the total mass distribution; hence, they might be concentrated 
toward the center of the galaxy cluster. In the same context. White (1987) suggested that 
cD envelopes and the dihuse light in Coma have the same origin: they are composed of stars 
tidally stripped from galaxies during Coma’s collapse. So GCs might have been stripped as 
well, creating an IGC population. 

Alternatively, if IGCs do exist in galaxy clusters, they may have been formed in situ 
in scenarios such as biased GC formation. West (1993) argued that the number density 
of GCs in biased formation scenarios is extremely sensitive to the presence and amplitude 
of a background held, such as galaxy cluster-sized and galaxy-sized perturbations in the 
primordial matter distribution. GCs would have been much more likely to form in developing 
protogalaxies than in complete isolation, but a small number of IGCs is expected. In these 
scenarios, the number density of IGCs would be very small, and also might follow the galaxy 
cluster mass distribution. 

From a galaxy formation point of view, as GCs are thought to be among the oldest 
objects in the Universe, they provide useful information about the galaxy formation and 
evolution processes. By studying a globular cluster system (GCS) associated with a single 
galaxy, the history of the host galaxy is probed. The hypothesis of the existence of an IGC 
population has been used to gain insight into the high-S'Ar problem. West et ah (1995) 
argued that a high value of Fat in galaxies located near the centers of galaxy clusters is the 
result of the accretion of a number of IGCs. They assumed that an IGC population exists in 
all galaxy clusters, and that the number density of IGCs is concentrated toward the center 
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of the galaxy cluster. 

Recently, Cote et al. (2001) did a dynamical analysis of the GCS associated with M87. 
They concluded that the metal-rich GCs of M87 formed with the galaxy, and that the metal- 
poor GCs are gradually accreting onto M87 along Virgo’s principal axis. This conld therefore 
be an example of how the accretion of IGCs would influence the galaxy’s evolutionary history. 

Although some authors consider that the accretion of IGCs is capable of modifying 
the properties of the GCS of a galaxy, there are a number of others that defend the view 
that the greater part of a GCS is formed in sitn from the protogalactic clond (Blakeslee et 
al. 1997; Forbes, Brodie, & Grillmair 1997; Harris et al. 1998; McLanghlin 1999). So 
the possible influence of the accretion of IGCs in galaxy formation and evolutionary models 
is still an open issne. In this context, the existence or not of IGCs and, if they do exist, 
their abnndance and distribution become crucial ingredients in the development of a unihed 
galaxy evolution scenario. 

Even though the existence of IGCs is often assumed, there is little evidence to date to 
prove this. Do GCs form only during starbursts in galaxies or can they form in isolation? If 
IGCs exist in all galaxy clusters, their number and spatial distribntion are fundamental pieces 
of information for testing the different GC formation models. Furthermore, the influence of 
IGGs in the evolntion of elliptical galaxies is still not clear. Is accretion of IGGs important 
dnring the galaxy’s formation and evolutionary history? If this is the case, is it important 
for all elliptical galaxies, or only for those located in some privileged locations inside galaxy 
clusters? In this paper the existence of an IGG population in the Goma galaxy cluster has 
been tested and qnantihed using the snrface-brightness fluctuations (SBF) techniqne. 

Our purpose is to obtain information about the IGG population from the study of the 
unresolved point sources by means of SBF analysis. The measured SBF signal is produced 
mainly by faint galaxies and the IGG population (if snch a thing exists). For this reason, 
dnring SBF measurements, a good characterization of the faint end of the backgronnd dif¬ 
ferential galaxy nnmber connts, n{m), is fnndamental. A nnmber of authors have measured 
the slope of the faint end of n{m) in the R hlter. Resnlts range from 7 = 0.39 (Tyson 1988) 
to 7 = 0.31-34 (Steidel & Hamilton 1993). This slope has also been measured in the Hubble 
Deep Field (Williams et al. 1996), with 7 = 0.36 for the V filter and 7 = 0.31 for the / hlter 
down to magnitnde 26. For the intermediate R hlter, the slope must be between 7 = 0.31 
and 0.36. From the SBF analysis presented in this paper, we provide not only information 
about the existence or not of an IGG popnlation, but also an estimate of the slope of the 
faint end of n{m). 
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2. Observations and data reduction 

Observations of the 16 giant elliptical galaxies and the 4 “blank” regions studied in this 
paper were done on 2000 April 25 and 27, with the 2.5 m Isaac Newton Telescope (INT) 
at the Roque de los Muchachos Observatory (La Palma), using the Wide Field Camera and 
the Sloan R Filter. This work is based on the results obtained by Marin-Franch & Apari- 
cio ( 2002 ), where the observations, photometric calibration and data reduction procedures 
are described in detail, as well as the SBF analysis. We provide here only a short outline of 
the latter. 

The concept of SBF was introduced by Tonry & Schneider (1988), who noted that, 
in the surface photometry of a galaxy too far to be resolved, a pixel-to-pixel fluctuation is 
observed due to the Poisson statistics of the spatial distribution of stars, GCs, background 
galaxies, etc. The variance of the fluctuation depends on the stellar population, the GCS, 
background galaxies, foreground stars, and, of course, the distance. It can be assumed that 
the total pixel-to-pixel variance of an image is the sum of all the independent contributions. 

The SBF technique involves spectral analysis of the signal and provides as a result the 
total PSF-convolved variance (Pq), produced by all objects whose spatial flux distribution is 
convolved with the PSF: 

-Po = o-gp + o-QC + ^BG + (1) 

where cr^p, ctqq, (t|q, and are the variances produced by the stellar populations, GCs, 
background galaxies, and foreground stars, respectively. If a population of IGCs is being 
analyzed, a^p can be neglected, and hence 

Pq = O'lGC + ^BG + *^1) (2) 

where (Tiqq is the pixel-to-pixel variance produced by the IGC population where this exists. 

The pixel-to-pixel variance produced by a class of object can be evaluated as the second 
moment of the differential number counts of that object population. In this way, (t|q can be 
obtained from the background n{m), htted to the resolved point sources, and (j| from the 
predicted star counts, computed making use of a Galaxy model. 

Once and are estimated and Pq is measured, can be easily obtained from 
equation 2 . Note that in Coma, GCs cannot be directly resolved from ground-based ob¬ 
servations (except the brightest individuals), so, in order to determine the surface number 
density of IGCs, SBF analysis becomes necessary. 

The surface number density of IGCs {Niqq) can be estimated from assuming a 
hxed GC luminosity function ( 99 ). A hxed 99 provides a hxed (Tjqq/Aigc ratio that can be 
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used to transform variances into surface densities. This ratio depends on the photometric 
calibration. 

Marin-Franch & Aparicio’s (2002) results used here are listed in Table 2. The name 
of the region studied is listed in column 1. In column 2, the distance to the Coma galaxy 
cluster center is also shown. Finally, the obtained Pq in the outer region of each galaxy (the 
region further away from the galactic center, where galactic GCs no longer exist) and in the 
“blank” areas are listed in column 3. The following discussion is based on these results. 


3. Results 

Let us consider regions close to the studied galaxies, but beyond the galactic halos. If 
IGCs existed then ctjqq would be 

<^1^00 = Po — (^bg ~ ^fs- (3) 

Pq is measured from SBF analysis. But, in order to obtain Niqq, <^bg <^1 must be 
deduced hrst. 


3.1. Contribution of foreground stars 

Coma is in a nearly perpendicular direction to the Galactic plane, so the expected 
number of foreground stars is very low and its contribution to the background variance is 
negligible when compared with the variance produced by the GCSs of the Coma galaxies 
(Marin-Franch & Aparicio 2002). But foreground stars can produce a signihcant contribu¬ 
tion to the total variance in intergalactic regions, so they must be considered if the variance 
produced by IGCs and the faint end of n{m) is to be analyzed. 

Bahcall & Soneira (1981) computed the predicted star counts in selected helds and 
photometric bands using their model for the Galaxy. In particular, they computed the 
predicted star counts in the northern Galactic pole in the R hlter down to magnitude 26.5 
(Fig. 1). By computing the second moment of these star counts, can be obtained for the 
North Galactic Pole. The result is 

.| = (l,50±0.15)xl0-(^)^ (4) 

assuming a 10% uncertainty for the star counts model predictions. As Coma is very close to 
the North Galactic Pole, we adopt this value and assume that it is constant over Coma. 
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3.2. Faint end of the differential galaxy number counts, n{m) 

The background galaxy variance, has to be obtained from the differential galaxy 
number counts, n{m). Tyson (1988) obtained a slope of 7 = 0.39 for the faint end of n{m) 
in filter R. More recently, Steidel & Hamilton (1993) obtained slopes in the range 7 = 
0.30-0.34, also in R. Unfortunately, the adopted value of 7 has considerable effects in the 
resulting IVigc- For this reason, we have evaluated iVicc considering three values for the 
n{m) slope, 7 = 0.30, 7 = 0.34, and 7 = 0.39, and has analyzed the results separately. The 
corresponding zero points of n{m) have been obtained from fits to the resolved point sources 
measured by Marin-Franch & Aparicio (2002) in which the slope is a fixed parameter. 

As an example, n{m) values are shown in Fig. 2 for the three slopes 7 = 0.30, 7 = 0.34, 
and 7 = 0.39 overplotted on the resolved point sources in region “Blank 3”. It can be seen 
that the three values of 7 are in principle compatible with the resolved galaxy distribution 
in this region. 

The variance (t|q can be now deduced as the second moment of n(m), and (Tiqc can be 
obtained from equation 3. In Table 2, and for each value of 7 , the results associated with each 
galaxy for (t|q are listed in column 4. The error bars of have been computed taking into 
account only the n{m) fitting uncertainty. Note that in Marin-Franch & Aparicio (2002) the 
error bar of was determined considering possible variations in the slope of the faint end 
of the n{m) (from 7 = 0.30 to 7 = 0.40), but in the present study, these possible variations 
are being considered explicitly by studying the cases 7 = 0.30, 0.34, and 0.39 separately. 


3.3. Intergalactic globular clusters 

Equation 3 can now be used to derive <Jiqq for each region. The results are given in 
column 5 of Table 2 . The number population of IGC, A^igC) can be estimated from Ujqq 
using the ctjqq/Aigc ratios listed column 6 . These ratios were computed in Marin-Franch 
& Aparicio (2002). The resulting A^igc is listed in column 7. Note that A^igg is expresed 
in IGCs/{"Y, while the rest of columns are expresed in pixels. The (Tiqq uncertainties 
have been computed adding in quadrature the uncertainties of Pq and (j|q. On the other 
hand, the uncertainties of the ratio ctjqc/A^igg have been neglected. The reason is that p 
is assumed to be universal (we assume the same cp parameters adopted in Marin-Franch 
& Aparicio 2002, a = 1.40 and = 27.42). In this sense, the uncertainty in the p 
parameters would introduce a rescaling factor common to all the results, but the relative 
errors would not change and the proofs of IGCs detection would remain valid. The origin 
of the galaxy-to-galaxy variations in the ratio (Tjqc/A^igg relies on the differences of the 
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photometric calibration constants. 

The values obtained for Niqq for the held regions associated with the 16 galaxies studied 
in Marin-Franch & Aparicio (2002) (dots) and for four additional intergalactic held regions 
(stars) are plotted in Fig. 3 versus the distance R to the Coma galaxy cluster center. Each 
panel corresponds to a value of 7 . Regions associated with galaxies were selected in Marin- 
Franch & Aparicio (2002) with the aim of sampling the background. They are beyond the 
GCS of each galaxy and they consequently sample the intergalactic held population. The 
four intergalactic “blank” regions are far from any galaxy. The fact that the average of Ajgc 
for these regions is the same as for the regions associated with galaxies conhrms that the 
held is really being sampled in both cases. 


The most signihcant characteristic of the plots in Fig. 3 is the hat distributions of values. 
The only point perhaps departing from this trend is that associated with the central galaxy, 
NGC 4874. Furthermore, it can be seen that Aigc is negative in almost all cases for 7 = 0.39, 
indicating that this value is unlikely. On the other hand, if 7 = 0.30, then A^igc is clearly 
a positive number, implying that an IGG population should exist. This population is also 
detected if 7 = 0.34, but less clearly. It should be noted that, if 7 = 0.30-0.34, Aigg does 
not correlate with R, and hence with the environment, except for the fact that the density 
associated with the central galaxy is higher. Rather, the IGG population would be uniformly 
spread all over Goma. We will discuss this issue further in Section 4. 


For each value of 7 , we have computed the average (Aigc) in Goma using all but the 
central points. The results are listed in Table 3 and have been plotted in Fig. 3 by dashed 
lines. Dotted lines show the 1 Upoints intervals (cipoints is the root mean square of the data 
point distribution). The estimated errors of the mean have also been computed using the 
usual formula 




^points 

v^A^’ 


( 5 ) 


and are the errors quoted in Table 3. 


4. Discussion 

4.1. The faint end of the differential galaxy number counts, n{m) 

The extreme cases considered above for the n{m) slope, 7 = 0.30 and 7 = 0.39, charac¬ 
terize two distinct scenarios. We have seen that the latter produces negative Aigc all over 
Goma, which indicates that n{m) must be shallower. 

Let us consider now in more detail the possibility 7 = 0.30. In this case. Fig. 3 shows 
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that IGCs would exist and would be spread all over Coma. Furthermore, no relation with 
the distance to the galaxy cluster center is apparent. However, all the GC formation sce¬ 
narios mentioned in §1 imply that Niqq should follow the galaxy cluster mass distribution. 
Assuming this and in order to test our results, we will consider three different galaxy clus¬ 
ter simplistic mass distributions, p{R), and analyze the expected A^igc distribution for each 
model. 


Makino (1994) has given for Coma a cut-off radius = 3.5°, a core radius = 0.24° 
(i.e. rt/r'c = 14.2) and, using iJo = 72 ± 8 km s“^ Mpc“^ (Freedman et ah 2001) for 
the Hubble parameter, a total mass 1.4 x lO^^M©. This value for the Hubble parameter 
corresponds to a distance of 97.2 Mpc and to = 5.9 Mpc and Tc = 0.4 Mpc. In this 
context, the three considered galaxy cluster mass distribution models correspond to a King 
model, a homogeneous sphere of radius equal to the cut-off radius of Coma {rt = 210') and 
a homogeneous sphere of radius equal to twice the cut-off radius of Coma. All them are 
plotted in the upper panel of Figure 4. All the mass distributions are normalized to a total 
galaxy cluster mass of 1.4 x IO^^Mq (Makino 1994). The vertical solid line represents the 
Coma cut-off radius, r*. 


Once a mass distribution model is assumed, the expected A^igc can be deduced pro¬ 
jecting the galaxy cluster spatial mass distribution and taking into account a GC formation 
efficiency. The latter is a measure of the fraction of mass converted into GCs. McLaugh¬ 
lin (1999) dehned the GC formation efficiency in galaxies as 

PGC 

cgc —-—- 

Pgas “r Pstars 

where pgC; Pgas, and Pstars are the mass densities of GCs, gas, and stars in the galaxy, 
respectively. McLaughlin (1999) found this definition of ecc to be universal from galaxy to 
galaxy and equal to 0.0026 ± 0.0005. 

We will use here a similar dehnition for the intergalactic medium that, for convenience, 
we dehne in terms of surface density: 

SiGC 
eiGC - — 

^tot 


( 6 ) 


where Sigc is the surface mass density of IGCs and Stot is the baryonic surface mass density 
of the galaxy cluster obtained projecting p{R). 


The relation between the baryonic surface mass density and the ICC surface number 
density is given by 

CGG^tot 

JVigG = - 


{m)Gc 


( 8 ) 
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where {m)GC is the mean mass of a GC. Assuming (m)Gc = 2.4 x 10^ (McLaughlin 1999), 
the lower panel of Fig. 4 shows A^igc/ggc as a function of R for the three mass distribution 
models plotted in the upper panel. This hgure shows that the baryonic mass distribution 
model required to produce an approximately constant, non-null, A^igc distribution along 
Coma (up to i? ~ 200', at least) should be similar to a very extended homogeneous sphere. 
The unrealistic nature of this model indicates that a value as low as 7 = 0.30 for the slope 
of the faint end of n{m) is unlikely. Moreover, it suggests that the IGC surface density, 
^iGC) is actually zero all over Coma, with the probable exception of the surroundings of the 
central galaxy. The value of 7 that reduces the mean (A^igc) to zero (excluding the central 
value) can be calculated from interpolation in the data listed in Table 3 and turns out to be 
7 = 0.36 ±0.01. 

The former is the slope of n{m) beyond the detection limit of i? = 23.5 reached by 
Tyson (1988) and Steidel & Hamilton (1993) and is valid at least down to the R magnitude 
signihcantly contributing to the SBF. This magnitude can be evaluated by considering that in 
Table 2 the error bars associated with are, on average, about 5%. If (t|q is computed for 
different limiting magnitudes, it is found that the variance produced by background galaxies 
fainter than rriR = 26.5 is lower than this value of 5%. The signal produced by fainter 
galaxies remains below the error bars. In conclusion, rriR = 26.5 can be assumed the faintest 
limit of our determination, so that the slope of the faint end of n{m) is 7 = 0.36 ±0.01 (this 
slope is valid down to mR = 26.5 at least). 

For comparison, it is interesting to note that Williams et al. (1996) have measured the 
slope of n{m) for deeper V and I magnitudes from the Hubble Deep Field. They obtain 
7 = 0.36 for V and 7 = 0.31 for I for the magnitude interval [23,26], and 7 = 0.17 for V and 
7 = 0.18 for I for the magnitude interval [26,29]. 


4.2. Intergalactic globular clusters 

In §4.1 we concluded that the most likely value for the general IGG density in Goma 
is zero, and that this is accomplished with a slope 7 = 0.36 ± 0.01 for the faint end of 
the background differential galaxy number counts, n{m). We will now further discuss the 
implications of these results on the IGG population and the overall baryonic mass distribution 
in Goma. 

In Fig. 5, Aigc obtained with 7 = 0.36 is plotted as a function of R, the distance to 
the Goma center. The points have been computed as an interpolation between the 7 = 0.34 
and 7 = 0.39 results listed in Table 2. It can be seen that the mean (Wgg) is zero, as 
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required. Only NGC 4874 seems to be surrounded by a significant number of IGCs. A King 
model, such as that plotted in Fig. 4, can now be compared with the IGG distribution. The 
efficiency eicc can be calibrated from the central IGG density Aigc = 0.059 and is 
Ggc = 7.01 X 10“®. Incidentally, this value is two orders of magnitude smaller than that 
obtained by McLaughlin (1999) for GGs belonging to galaxies. 

The former ejGc aiid the concentration parameter rt/r^ = 14.2 derived by Makino (1994) 
for the intergalactic matter in Goma associated with the X-ray radiation can be used to plot 
the expected distribution of IGGs. The result is shown in Fig. 5, upper panel (solid line). 
This would be the distribution of IGG in Goma if; (i) the X-ray emission properly samples 
the baryonic intergalactic mass distribution in Goma; (ii) the GG excess observed in the 
central region is caused by a genuine IGG population bound by the general Goma potential 
well, and (hi) if the efficiency e is constant all over Goma. This distribution lies beyond 
the error bars of several single Aigc values. A higher concentration parameter is required 
to properly fit the Aigc point distribution. The best fit to a King distribution produces 
r^/rc = 81.5 and eiGc = 9-5 x 10“^. The result has been plotted in the lower panel of Fig. 5. 

The reduced chi-square for the best fit results xt = 0-84 which indicates that the fit is 
good enough. On the other hand, xt = 3.0 for the rt/rc = 14.2 king model fit. For z/ = 19 
(as in our case), the integral probability of the xl distribution is 0.001 for xt > 2.31, showing 
that the rt/rc = 14.2 king model fit can be rejected with a high confidence level. 

The concentration parameter rt/rc = 81.5, implies a core radius of = 154". This is 
of the order of, or smaller than, the size of the Goma central galaxy NGG 4874. In fact, in 
Marin-Franch & Aparicio (2002), the spatial structure of the NGG 4874 GGS was studied 
out to a distance of 161.4" from the galaxy center, and its edge was probably not reached. 
This suggests that the obtained A^igg excess around the center of Goma is produced by the 
outermost extension of the NGG 4874 primitive GGS rather than being a trace of an IGG 
system. In other words, the spatial distribution of A^igg shown in Fig. 5 is probably the 
map of the NGG 4874 GGS. 

This result suggests that an IGG population does not exist in Goma. This points 
toward GGs having been formed only, or almost only, from protogalactic clouds. None of 
them, or perhaps very few, could have formed in isolated regions. This is in good agreement 
with the in situ formation scenario for the GGSs (Blakeslee et al. 1997; Forbes, Brodie, & 
Grillmair 1997; Harris et al. 1998; McLaughlin 1999). If the distribution of IGGs followed 
the distribution of dark matter in galaxy clusters (Muzzio 1987; West 1993), it would be 
reasonable to expect that Wgg would show a spatially extended gradient across de galaxy 
cluster. This not being the case, it seems inappropriate to support the relationship between 
IGG and dark matter distributions, although it is true that the relationship could still exist 
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but not be strong enough to be detected. 

On the other hand, a number of authors have argued that accretion of IGCs might infln- 
ence the formation and evolution of galaxies (Muzzio 1987; White 1987; West et ah 1995; 
Cote et al. 2001). Our results indicate that accretion of IGCs is not a signihcant effect in 
galaxy formation and evolutionary processes in the Coma galaxies, since our conclusion is 
that IGCs do not exist in Coma. 

Gnrzadyan & Mazure (2001) fonnd the existence of three snbgronps of galaxies in Coma, 
one of them associated with the cD galaxy NGC 4874 and the other two with NGC 4889 
and NGC 4839. They conclude that the non-stationarity of the dynamical processes at work 
in the Coma core is due to the merging of small-scale groups of galaxies. In this context, 
each subgronp formed separately and then the merger between the different gronps took 
place. In Marin-Franch & Aparicio (2002), hve galaxies belonging to the Gnrzadyan & 
Maznre (2001) snbgronp 1 (the NGC 4874 group) were analyzed showing that a relation 
between Sn and the distance to the central galaxy NGC 4874 was compatible with the data. 
The possibility of accretion of IGCs being responsible for this relation was discussed. In the 
light of our new results and owing to the absence of IGCs found here, this possibility should 
be also rejected. 


5. Conclusions 

The existence of an ICC population in the Coma galaxy cluster has been tested using the 
SBF. The measured SBF signal is mainly produced by faint galaxies and the ICC population 
(if it exists). For this reason, the Wgc results reached and the n{m) faint end are connected. 
As we do not have a precise enough measurement of n(m) in the Coma galaxy cluster, Wcc 
has been obtained considering three different valnes for the n{m) slope (7 = 0.30, 0.34, and 
0.39) and the resnlts have been analyzed separately. The main conclusions of this stndy are 
snmmarized here: 

• For the three 7 values studied, the results obtained for Wgc have been shown not to 
correlate with the distance to the central galaxy NGC 4874, and hence with the envi¬ 
ronment. Comparing the Wgc results obtained for each valne of 7 with the different 
Coma galaxy cluster mass distribution model predictions, we suggest that Aigc must 
be zero all over Coma, with the probably exception of the surronndings of the central 
galaxy. The most likely n{m) slope that makes (Aigc) = 0 is 7 = 0.36 ± 0 . 01 . This 
slope has been shown to be valid up to mu = 26.5. For comparison, previous results 
for the n{m) slope in the R hlter range from 7 = 0.39 (Tyson 1988) to 7 = 0.31-0.34 
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(Steidel & Hamilton 1993). This slope has also been measured in the Hubble Deep 
Field North in the magnitude interval [23,26], resulting in 7 = 0.36 for the V hlter and 
7 = 0.31 for the I hlter (Williams et ah 1996). Our new 7 determination falls between 
previous measurements. 

• Our results seem to be in good agreement with the in situ formation for the GCSs 
(Blakeslee et ah 1997; Forbes, Brodie, & Grillmair 1997; Harris et ah 1998; McLaugh¬ 
lin 1999). The former conclusions point toward GGs having been formed only, or 
almost only, from protogalactic clouds. None of them, or perhaps very few, could have 
formed in isolated regions. 

• Muzzio (1987) and West (1993) argued that there might exist a relationship between 
IGG and dark matter distributions. In the light of our results it seems inappropriate 
to support this idea, although it is true that the relationship could still exist but not 
be strong enough to be detected. 

• Finally, a number of authors have argued that accretion of IGGs might influence the 
formation and evolution of galaxies (Muzzio 1987; White 1987; West et al. 1995; Gote 
et al. 2001). Our results suggest that accretion of IGGs might not be a signihcant 
effect in galaxy formation and the evolutionary processes in the Goma galaxies since 
our conclusion is that IGGs do not exist in Goma. 


This article is based on observations made with the 2.5 m Isaac Newton Telescope 
operated on the island of La Palma by the ING in the Spanish Observatorio del Roque 
de Los Muchachos. This research has been supported by the Institute de Astrofisica de 
Ganarias (grant P3/94), the DGESIG of the Kingdom of Spain (grant PI97-1438-G02-01), 
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Fig. 2.— Background differential galaxy number counts, n{m), fitted to resolved point source 
distribution in the region “Blank 3”. The £t has been performed using the three values for 
the n{m) slope 7 = 0.30, 0.34, and 0.39 discussed in the text as fixed parameters. The 
error bars represent the Poissonian error of the detection counts. Solid line represents the 
obtained n{m) in this paper (see §4.1). 
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Fig. 3.— The IGC population in Coma plotted for three possible values of the n{m) slope, 
7 . Horizontal axis represents the distance to the Coma galaxy cluster center, and vertical 
axis represents the ICCs surface number density. Full lines indicate the zero density level. 
Dashed lines represent the mean for each case while dotted lines show the iln of the mean. 
Full circles represent the outskirts of the 16 galaxies studied, while open stars represent the 
four control fields. See text for details. 
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Fig. 4.— Upper panel: The total mass density (p) is plotted as a function of the distance 
to the galaxy cluster center (R) for three different galaxy cluster mass density distribu¬ 
tions; King (1962) model (long dash), homogeneous sphere of radius equal to the Coma 
radius (dots) and homogeneous sphere of radius twice the Coma radius (short dash). Lower 
panel: ICC surface number density normalized to the globular cluster formation efficiency 
(^igc/ugc) for each of the former mass distributions. See text for details. 
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Fig. 5.— iViGc interpolated for the resulting case 7 = 0.36. Upper panel; The theoretical 
expectation from the King (1962) model has been htted to our central Niqq result consid¬ 
ering the Makino (1994) core radius (solid line); Lower panel: iViGc results have been htted 
with a King (1962) model considering rt/r^ and cigc as free parameters. The best ht is 
plotted with a solid line. 
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Table 1. Notation summary. 


Globular Cluster GC 

Intergalactic Globular Cluster IGC 

Globular Cluster System GCS 

Surface Brightness Fluctuations SBF 

Differential Galaxy Number Counts n(m) 

Differential Galaxy Number Counts Slope 7 

Globular Cluster Luminosity Function (/? 





Table 2 . IGCs in the environments of the 16 Coma galaxies^?^ 


Galaxy 

R 

Po 


2 

°'lGC 


^IGC 


(') 

10“"^ X f 

\ sxpix J 

lO”"' X f = ) 

\ sxpix J 

lO”"' X f ) 

\ sxpix J 

^igc/^igc('’| 

IGCs/(")^ 

7 = 0.30 

NGC 4874 

0 

31.8 ± 1.2 

22.6 ± 1.2 

7.7 ± 1.7 

0.067 

0.104 ± 0.023 

NGC 4889 

8.05 

26.2 ± 2.1 

19.8 ± 1.4 

4.9 ± 2.9 

0.067 

0.066 ± 0.039 

IG 4012 

10.68 

18.0 ± 1.1 

14.9 ± 1.2 

1.6 ± 1.6 

0.067 

0.022 ± 0.022 

IG 4021 

10.88 

16.4 ± 1.0 

12.5 ± 0.8 

2.4 ± 1.3 

0.067 

0.032 ± 0.017 

IG 4026 

12.62 

17.2 ± 1.1 

12.0 ± 0.9 

3.7 ± 1.4 

0.067 

0.049 ± 0.018 

IG 4041 

16.38 

21.2 ± 1.1 

15.6 ± 1.1 

4.1 ± 1.6 

0.067 

0.055 ± 0.021 

IG 4045 

19.85 

19.7 ± 1.1 

15.3 ± 0.9 

2.9 ± 1.4 

0.067 

0.039 ± 0.019 

IG 4051 

19.95 

26.1 ± 2.1 

20.4 ± 0.8 

4.2 ± 2.2 

0.067 

0.056 ± 0.029 

IG 3976 

6.78 

11.5 ± 3.0 

8.4 ± 1.3 

1.6 ± 3.2 

0.045 

0.032 ± 0.064 

IG 3959 

12.67 

15.2 ± 1.0 

13.6 ± 0.8 

0.1 ± 1.3 

0.059 

0.002 ± 0.026 

MGG +5 -31 -063 

8.93 

13.5 ± 0.6 

12.3 ± 0.3 

-0.3 ± 0.7 

0.059 

-0.005± 0.012 

NGC 4839 

43.12 

10.3 ± 1.0 

8.7 ± 0.1 

0.1 ± 1.0 

0.046 

0.002 ± 0.020 

NGC 4840 

37.43 

9.5 ± 0.8 

8.6 ± 0.3 

-0.6 ± 0.9 

0.046 

-0.012± 0.018 

NGC 4816 

52.60 

16.5 ± 1.1 

10.8 ± 0.4 

4.2 ± 1.2 

0.060 

0.064 ± 0.018 

NGC 4673 

217.07 

15.5 ± 0.6 

12.9 ± 0.4 

1.1 ± 0.7 

0.067 

0.015 ± 0.009 

IG 3651 

290.32 

16.7 ± 0.6 

12.9 ± 0.3 

2.3 ± 0.7 

0.067 

0.031 ± 0.009 

Blank 1 

40 

10.3 ± 0.3 

6.7 ± 0.4 

2.1 ± 0.5 

0.060 

0.031 ± 0.007 

Blank 2 

56 

14.5 ± 2.0 

9.8 ± 0.3 

3.2 ± 2.0 

0.060 

0.048 ± 0.030 

Blank 3 

80 

11.0 ± 2.0 

7.7 ± 0.3 

1.8 ± 2.0 

0.060 

0.027 ± 0.030 

Blank 4 

280 

16.0 ± 3.0 

11.4 ± 0.5 

3.1 ± 3.0 

0.060 

0.046 ± 0.044 


7 = 0.34 


NGC 4874 

0 

31.8 ± 1.2 

24.7 ± 0.8 

5.6 ± 1.4 

0.067 

0.075 ± 0.019 

NGC 4889 

8.05 

26.2 ± 2.1 

21.8 ± 1.2 

2.9 ± 2.4 

0.067 

0.039 ± 0.032 

IG 4012 

10.68 

18.0 ± 1.1 

15.7 ± 1.2 

0.8 ± 1.6 

0.067 

0.011 ± 0.022 

IG 4021 

10.88 

16.4 ± 1.0 

13.5 ± 0.9 

1.4 ± 1.3 

0.067 

0.019 ± 0.018 

IG 4026 

12.62 

17.2 ± 1.1 

14.4 ± 0.9 

1.3 ± 1.4 

0.067 

0.017 ± 0.018 

IG 4041 

16.38 

21.2 ± 1.1 

17.1 ± 1.0 

2.6 ± 1.5 

0.067 

0.035 ± 0.020 

IG 4045 

19.85 

19.7 ± 1.1 

16.8 ± 1.0 

1.4 ± 1.5 

0.067 

0.019 ± 0.020 

IG 4051 

19.95 

26.1 ± 2.1 

22.3 ± 0.8 

2.3 ± 2.2 

0.067 

0.031 ± 0.029 

IG 3976 

6.78 

11.5 ± 3.0 

9.2 ± 1.3 

0.8 ± 3.2 

0.045 

0.016 ± 0.064 

IG 3959 

12.67 

15.2 ± 1.0 

15.0 ± 0.8 

-1.3 ± 1.3 

0.059 

-0.020 ± 0.020 

MGG +5 -31 -063 

8.93 

13.5 ± 0.6 

13.4 ± 0.3 

-1.5 ± 0.7 

0.059 

-0.023 ± 0.011 

NGC 4839 

43.12 

10.3 ± 1.0 

9.6 ± 0.1 

-0.8 ± 1.0 

0.046 

-0.016 ± 0.020 

NGC 4840 

37.43 

9.5 ± 0.8 

9.4 ± 0.3 

-1.4 ± 0.9 

0.046 

-0.027 ± 0.017 

NGC 4816 

52.60 

16.5 ± 1.1 

13.8 ± 0.3 

1.2 ± 1.1 

0.060 

0.018 ± 0.016 

NGC 4673 

217.07 

15.5 ± 0.6 

13.7 ± 0.3 

0.3 ± 0.7 

0.067 

0.004 ± 0.009 

IG 3651 

290.32 

16.7 ± 0.6 

14.1 ± 0.2 

1.1 ± 0.6 

0.067 

0.015 ± 0.008 

Blank 1 

40 

10.3 ± 0.3 

7.4 ± 0.4 

1.4 ± 0.5 

0.060 

0.021 ± 0.007 

Blank 2 

56 

14.5 ± 2.0 

12.2 ± 0.3 

0.8 ± 2.0 

0.060 

0.012 ± 0.030 

Blank 3 

80 

11.0 ± 2.0 

8.7 ± 0.3 

0.8 ± 2.0 

0.060 

0.012 ± 0.030 

Blank 4 

280 

16.0 ± 3.0 

12.5 ± 0.5 

1.0 ± 3.0 

0.060 

0.015 ± 0.045 


7 = 0.39 


NGC 4874 

0 

31.8 ± 1.2 

27.7 ± 0.6 

2.6 ± 1.3 

0.067 

0.034 ± 0.017 

NGC 4889 

8.05 

26.2 ± 2.1 

24.3 ± 1.2 

0.4 ± 2.4 

0.067 

0.005 ± 0.030 

IG 4012 

10.68 

18.0 ± 1.1 

16.9 ± 1.2 

-0.4 ± 1.6 

0.067 

-0.005 ± 0.020 

IG 4021 

10.88 

16.4 ± 1.0 

14.9 ± 1.1 

0.0 ± 1.5 

0.067 

0.000 ± 0.019 

IG 4026 

12.62 

17.2 ± 1.1 

17.8 ± 0.9 

-2.1 ± 1.4 

0.067 

-0.028 ± 0.019 
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Table 2—Continued 


Galaxy 

R 

(') 

Po 

^BG 

10-4 .. f 

2 


Nigg 

IGCs/(")2 

'' sxpix J 

'' sxpix J 

' ' sxpix J 

IC 4041 

16.38 

21.2 ± 1.1 

19.2 ± 1.0 

0.5 ± 1.5 

0.067 

0.007 ± 0.021 

IC 4045 

19.85 

19.7 ± 1.1 

18.6 ± 1.2 

-0.4 ± 1.6 

0.067 

-0.005 ± 0.020 

IC 4051 

19.95 

26.1 ± 2.1 

24.9 ± 0.9 

-0.3 ± 2.3 

0.067 

-0.004 ± 0.030 

IC 3976 

6.78 

11.5 ± 3.0 

10.3 ± 1.4 

-0.4 ± 3.3 

0.045 

-0.008 ± 0.066 

IC 3959 

12.67 

15.2 ± 1.0 

16.7 ± 0.7 

-3.0 ± 1.2 

0.059 

-0.046 ± 0.018 

MCG +5 -31 -063 

8.93 

13.5 ± 0.6 

12.1 ± 0.3 

-0.1 ± 0.7 

0.059 

-0.002 ± 0.014 

NGC 4839 

43.12 

10.3 ± 1.0 

10.7 ± 0.1 

-1.9 ± 1.0 

0.046 

-0.037 ± 0.019 

NGC 4840 

37.43 

9.5 ± 0.8 

10.5 ± 0.2 

-2.5 ± 0.8 

0.046 

-0.049 ± 0.016 

NGC 4816 

52.60 

16.5 ± 1.1 

18.5 ± 0.6 

-3.5 ± 1.3 

0.060 

-0.052 ± 0.019 

NGC 4673 

217.07 

15.5 ± 0.6 

14.8 ± 0.2 

-0.8 ± 0.6 

0.067 

-0.011 ± 0.008 

IC 3651 

290.32 

16.7 ± 0.6 

15.7 ± 0.2 

-0.5 ± 0.6 

0.067 

-0.007 ± 0.008 

Blank 1 

40 

10.3 ± 0.3 

8.5 ± 0.4 

0.3 ± 0.5 

0.060 

0.004 ± 0.006 

Blank 2 

56 

14.5 ± 2.0 

15.2 ± 0.3 

-2.2 ± 2.0 

0.060 

-0.033 ± 0.030 

Blank 3 

80 

11.0 ± 2.0 

10.1 ± 0.3 

-0.6 ± 2.0 

0.060 

-0.009 ± 0.030 

Blank 4 

280 

16.0 ± 3.0 

14.1 ± 0.5 

0.4 ± 3.0 

0.060 

0.006 ± 0.045 


(“^Columns Pq, CgQ and (Tjqq/TVigc have been extracted from the analysis done in Marin-Franch & Aparicio (2002). 1 

pix = 0.333^^ 

Theoretical value based on the assumption that the p is universal with cr = 1.40 and = 27.42. This value depends on 
the photometric calibration. 
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Table 3. Mean surface number density of IGCs in Coma. 


7 

(Mgc) 


(iGCs/r)b 

0.30 

0.032 ± 0.005 

0.34 

0.010 ± 0.004 

0.39 

-0.014 ± 0.004 



